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We investigate the range of validity of the recently developed steady-state ab-initio laser theory 
(SALT). While very efficient in describing various microlasers, SALT is conventionally believed not 
to be applicable to lasers featuring fully or nearly degenerate pairs of resonator modes above the 
lasing threshold. Here we demonstrate how SALT can indeed be extended to describe such cases as 
well, with the effect that we significantly broaden the theory’s scope. In particular, we show how to 
use SALT in conjunction with a linear stability analysis to obtain stable single-mode lasing solutions 
that involve a degenerate mode pair. Our flexible and efficient approach is tested on one-dimensional 
ring lasers as well as on two-dimensional microdisk lasers with broken symmetry. 
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INTRODUCTION 

Microcavity lasers are essential elements in modern 
photonics and have been realized with cavities of very 
different shape and with various lasing mechanisms [1- 
12]. The nonlinear lasing behavior of these systems 
can, in principle, be modeled using the semi-classical 
Maxwell-Bloch (MB) equations [13-16]. Due to their 
time-dependent nature these equations are, however, usu¬ 
ally difficult to solve for all but the most simple cases. 
In recent years, a much more efficient approach named 
steady state ab initio lasing theory (SALT) has emerged, 
which can be used to describe the steady-state las¬ 
ing of lasers [17-23]. Among other advances, this new 
framework has shed light on weakly-scattering random 
lasers [24], on pump-induced exceptional points [11, 12, 
25, 26] and on coherent perfect absorption [27, 28] and 
has opened up new ways of controlling the emission pat¬ 
terns of random as well as of microcavity lasers [29, 30]. 
One of the major drawbacks of SALT is that its conven¬ 
tional formulation fails for the simulation of microlasers 
with nearly degenerate modes as occurring, e.g., in whis¬ 
pering gallery mode resonators with an inherent symme¬ 
try [2, 6, 8, 10, 12, 31]. 

Here, we present an approach to generalize SALT to 
such cases. This extension allows us to observe, among 
other phenomena, that nearly degenerate modes may 
merge into a single-mode. These steady-state solutions 
are, however, not necessarily stable with respect to small 
time-dependent spatial perturbations. In particular, it 
has already been shown, that in highly symmetric sys¬ 
tems such as ring lasers not every steady state solu¬ 
tion of the MB equations is necessarily stable [32, 33]. 
Hence, the stability of the solutions obtained from our 
extended SALT approach has to be verified. Up to now 
such additional stability checks were always done us¬ 
ing direct time-dependent simulations of the MB equa¬ 
tions [18, 25, 34, 35]. One of the reasons for using SALT, 
however, is exactly to avoid this kind of computationally 


very demanding numerics. 

In this work, we thus introduce a much more efficient 
way to determine the stability of the SALT solutions 
based on a rigorous linear stability analysis for single¬ 
mode steady-state solutions. Furthermore, our work ex¬ 
tends the scope of SALT to previously inaccessible pa¬ 
rameter regimes, where even bifurcating solutions (sta¬ 
ble or unstable) of the nonlinear equations can now be 
appropriately dealt with. 

SHORT REVIEW OF SALT 

In semi-classical laser theory, the dynamics of a laser 
is governed by the interaction of classical fields with an 
ensemble of two-level atoms as described by the so-called 
Maxwell-Bloch (MB) equations. Restricting the fields to 
ID or to the transverse magnetic (TM) polarization in 
2D, the electric field and polarization become scalar [36]. 
Using the rotating wave approximation, the MB equa¬ 
tions can be derived as follows [16] 

eP + = V 2 P + - P + 

P+ = - (lLJ a + 7 ± ) P+ - i^E+D 

£ t= D a -D + i{E^)’-c. c.). 

In these non-linear partial differential equations, P + (x, t) 
and P + (x,£) denote the positive frequency components 
of the electrical field and the polarization of the medium, 
respectively. The quantity D(x,t) is the population in¬ 
version of the two-level atoms and Dq(x) stands for the 
externally imposed pump strength. The constants in the 
MB equations describe the properties of the cavity and of 
the gain material: the dielectric function e(x), the tran¬ 
sition frequency of the two level atoms co a , as well as the 
decay rates of the polarization j± and of the population 
inversion qy. The boundary conditions for the equations 
above are typically outgoing boundary conditions, which 


2 


numerically can, e.g., be implemented using a perfectly 
matched layer [22, 37]. 

The natural units in Eqs. (1) and all example systems 
in this work can be converted back to SI units by choosing 
an appropriate length scale L, multiplying all lengths in 
the example systems by this quantity L, multiplying the 
quantities (J^, 7|| , y_L, w a ) by c/Z, where c is the speed of 
light, dividing ^ by L and finally by transforming E, P, 
and D as follows Et = P+ = P+ . _ , 

bI ?V7||7x 5 ^0-^7117-L ’ 
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and Psi = P ^ £q • The variable g is the transition 
dipole moment of the two-level atoms. 

For microlaser systems lasing in steady-state, only a fi¬ 
nite number of modes in the system are active. This is the 
case described by SALT, in which the MB equations are 
simplified to a set of time-independent, non-Hermitian, 
nonlinear and coupled Helmholtz equations. The solu¬ 
tions of these SALT equations are much more efficiently 
calculated than those of the MB equations [18, 25, 34, 35]. 

The cornerstone of the SALT equations is a multi- 
periodic ansatz for the electro-magnetic field as well as 
for the polarization 


pump strengths Do, the traditional SALT algorithm will 
not be applicable when both modes of such a pair move 
across the lasing threshold. For the completely degener¬ 
ate case this problem is even more acute. 

In the following, we will provide an extension of the 
SALT algorithm, which is able to overcome this signifi¬ 
cant drawback. It does so by taking into account that 
degenerate modes with a fixed relative phase can be ex¬ 
pressed as a single active lasing mode within the SALT 
formalism. For nearly degenerate modes, we find that 
such a mode pair can become dynamically stable in the 
form of a single lasing mode, allowing us to treat the 
solution again with the above SALT ansatz. In order 
to present our approach as clearly as possible, we will 
focus here on the case of single-mode lasing only, keep¬ 
ing in mind that our algorithm can be generalized to 
multiple pairs of degenerate modes or multiple modes in 
general [39]. 

For a single lasing mode the ansatz (2) satisfies the 
MB Eqs. (1) exactly, leading to the following single-mode 
SALT equation for the electric field E\ and for the real 
lasing frequency ay: 


N 

E+(x,t) = Y / M%) e“ iuvt 

M=1 ( 2 ) 
p+( X ,t) = j2w *-*”»*, 

n=l 

where each triplet (F^, P M , ay) represents a mode E M of 
the system lasing at the real frequency ay. Inserting the 
ansatz (2) into the last equation of (1) results in 

D = 7 || (D 0 -D)+ l -f - c.c.) (3) 

with terms of the form P u E^ l e l ^ UJu ~ UJ ^ t , which explicitly 
depend on time. For multiple lasing modes, D will thus 
never be completely static. However, if the timescale 
1 /(ay, —l o^) on which these terms oscillate is much shorter 
than the timescale on which D varies (1/yy), their con¬ 
tribution can be neglected. In other words, for systems 
where |c— oy| yy holds for all pairs fi, v of active 
modes [38] the inversion can be approximated to be sta¬ 
tionary [17, 22]. 

This stationary inversion approximation (SIA) is, how¬ 
ever, not well satisfied in macroscopic lasers with a large 
density of modes as well as in microcavity lasers with an 
inherent symmetry. We will focus here on the latter case 
and shall consider ring- or microdisk lasers with degener¬ 
ate modes or slightly perturbed versions of these systems 
with nearly degenerate modes. The frequency splitting of 
these nearly degenerate modes will typically violate the 
condition Ac o yy for realistic values of yy. Since not 
only the lasing frequencies of these modes are very close 
to each other but also their thresholds are at comparable 
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E 1 (x)=0, (4) 


where 


j- | . 5 

CJi - L0 a + zy_L 

and the boundary conditions are the same as for the 
MB equations. Note, that Eq. (4) depends nonlinearly 
both on the frequency ay as well as on the shape of the 
mode E\(x) through a self-saturation spatial hole burn¬ 
ing interaction. It can be straightforwardly solved using 
a Newton-Raphson solver as described in detail in [22]. 
The solver requires an initial guess which can be obtained 
by tracking the modes in the system from an initial value 
at zero pump strength up to the pump strength of in¬ 
terest. This procedure has the advantage that for small 
pump values (below the lasing threshold) one can use the 
following eigenvalue problem that is linear with respect 
to the mode profiles, 

V 2 + (e^ + fiDoOr))^ 2 

where T* = y±/(dy — uo a + iy_i_). (Note that we label 
all quantities below the lasing threshold by overbars.) 
For pump strength Do = 0 the complex eigenvalues ay 
have a negative imaginary part. When increasing the 
pump strength, the eigenfrequencies will typically move 
towards the real axis (interesting exceptions to this rule 
are discussed in [11, 12, 25]). Once the first eigenvalue ay 
crosses the real axis (we assume its index i is 1) it can be 
used as a guess for solving Eq. (4) for the first lasing mode 


Ei(x)= 0 , ( 6 ) 
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Ei with real frequency uj\. This solution can then be 
tracked to even higher pump strengths beyond the lasing 
threshold by repeatedly solving Eq. (4) for increasing Dq 
while using the solution from the previous step as an 
initial guess. Alternative methods to solve Eq. (4) based 
on an expansion of laser modes in a bi-orthogonal basis 
of “constant-flux states” [17] work in a similar way but 
will not be considered here. 

In the traditional SALT algorithm the validity of a 
single-mode solution {Ei,ui} of Eq. (4) is only indirectly 
assessed by keeping track of the remaining passive modes 
of the system. While only one mode {Ei, uji} is lasing all 
these other modes (with i ^ 1) have to solve the following 
nonlinear, non-Hermitian eigenvalue problem 


V 2 + e(x) + f i 


D 0 {x) 


i + ir^O)! 2 




Ei(x) = 0. (7) 


The single-mode solution is only valid as long as all other 
eigenmodes Ei have an eigenvalue ay with imaginary part 
less or equal to 0, i.e., 


Mi ^ 1 : Im(oy) < 0 (8) 


If, however, any other of these eigenvalues crosses the real 
axis, the corresponding eigenmode is assumed to be ac¬ 
tive and incorporated as an additional lasing mode into 
the active SALT equations [22]. As long as the pres¬ 
ence of this second lasing mode does not violate the SI A, 
the corresponding two-mode lasing solution is considered 
stable (as was previously verified using FDTD simula¬ 
tions [18, 34, 35]). 

As we will demonstrate through a comparison to time- 
dependent solutions of the MB equations, this simple cri¬ 
terion can not be applied for nearly degenerate modes. 
In particular, we show that stable single-mode solutions 
may exist even though one of the other eigenmodes in 
the system features an eigenvalue dy with positive imag¬ 
inary part. In order to be able to correctly determine 
the stability of a mode when using the SALT Eq. (4), we 
introduce below a rigorous stability criterion based on a 
linear stability analysis. 

One of the consequences of this new strategy is that 
eigenmodes have to be continuously tracked even when 
their eigenvalues cross the real axis without, however, in¬ 
cluding them as an active lasing solution. Doing this, 
we find that these modes exhibit complicated frequency 
shifts and bifurcations when varying the pump strength 
Dq. Examples of this kind will be discussed in the sub¬ 
sequent sections. 


EXAMPLE 1: SYMMETRIC ID RING LASER 

To convey an understanding of how a lasing system 
with degenerate modes can be described in the SALT 
framework, let us first consider the well-known example 


of a rotationally symmetric ring laser whose solution can 
become unstable in certain parameter regimes [32, 33, 
40-48]. We model the system as a one-dimensional, ho¬ 
mogeneous medium with periodic boundary conditions 
(see Fig. la for an illustration). To incorporate losses 
through absorption and outcoupling, we set the index of 
refraction to a complex value. 

Solving Eq. (6) for the unpumped system produces 
a set of eigenstates where a two-dimensional 

eigenspace is associated with every eigenvalue dy due to 
the rotational symmetry of the system. This eigenspace 
contains standing waves of the form e lUiX ±c.c. as well as 
traveling waves of the form e ±luJiX . While these two pairs 
of states as well as their superpositions solve Eq. (6) be¬ 
low threshold, the non-linear spatial hole-burning term in 
the SALT Eq. (4) prevents arbitrary superpositions from 
being valid solutions above the threshold. This leaves 
only two possible steady-state lasing solutions of the rin- 
glaser at dy, corresponding to the well known clockwise 
and counterclockwise traveling wave states. However, 
from the literature it is known that ring lasers show 
complex behavior, including the fact that these travel¬ 
ing wave solutions are not always stable [32, 33, 40-48]. 

Our goal here will be to find the single-mode solutions 
with SALT and to identify the corresponding regions of 
stability. To approach this problem first in the most 
general way (i.e., independently of the employed SALT 
approach), we set up a finite-difference time-domain 
(FDTD) method based on a Yee lattice [49, 50]. This 
tool allows us to solve the MB equations Eq. (1) directly, 
including the full temporal evolution starting from an ini¬ 
tial distribution of the electric field. Using this approach 
we first confirmed that in the single-mode regime indeed 
only the traveling modes are stable in certain parameter 
regimes. To assess the latter, the system was initialized 
in the traveling-wave solution obtained from the SALT 
Eq. (4) [51] and then left to evolve for a certain amount 
of time. We analyzed whether at later times the system 
remained in the same steady-state solution as at the be¬ 
ginning of the simulation. If the system stayed in the 
same state (and did not show any signatures of deviating 
from this state), we considered the solution to be stable. 

The stability diagram resulting from these FDTD sim¬ 
ulations is shown in Fig. lc. We find that the stabil¬ 
ity of the single-mode SALT solutions not only depends 
on the external pump strength D 0 , but also on the in¬ 
version decay rate yy. This finding concurs with pre¬ 
vious results [32] and shows that while the single-mode 
SALT solution is an exact solution of the MB equations, 
it is not necessarily a stable one. For systems with de¬ 
generate passive modes (as well as with closely spaced 
modes |Acj| < yy discussed below), the stability of any 
solution obtained from SALT is therefore not guaranteed 
and needs to be independently verified. Using FDTD for 
such a verification, as we did above, is, however, much 
too costly from a numerical point of view, in particu- 
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Figure 1. (Color online) Results for a one-dimensional ring 
laser system with circumference L = 1, 7_i_ = 1 and uj a — 61. 
The openness of the system is modeled by a lossy dielectric 
function yje(x) — 1 + 0.0002L (a) Sketch of the ring laser 
system, which can accommodate clockwise (CW) or counter¬ 
clockwise (CCW) traveling modes, (b) Lasing mode inten¬ 
sity versus applied pump strength for the traveling wave so¬ 
lution. Note that the gain parameters have been chosen such 
that only a single pair of degenerate modes reaches the lasing 
threshold. The traveling wave solution in which the system 
lases has a frequency of uj ~ 62.65 which stays constant while 
the pump strength is increased, (c) Stability analysis of the 
SALT results. Shown are the stability of this traveling wave 
solution under variation of the pump strength Do and the re¬ 
laxation rate of the inversion yy. The results of the FDTD 
simulation are color-coded with green (light gray) marking 
parameter combinations where the SALT solution is stable, 
whereas the red (dark gray) region represents unstable be¬ 
havior. Blue dashed lines show the independent results of the 
linear stability analysis, which provides an excellent estimate 
for the border between stable and unstable regions. 


lar, as this would nullify the computational advantages of 
SALT. Also the stability analysis presented in [32] for ro- 
tationally symmetric ring lasers is rather limited in scope, 
as it relies on the availability of exact analytic solutions 
for the passive modes. We will thus develop below a more 
general and rigorous framework to analyze the stability 
of SALT solutions that should be generally applicable. 


LINEAR STABILITY ANALYSIS 

This section contains an overview of a linear stability 
analysis for solutions of SALT Eq. (4) (a full derivation 
can be found in appendix A). Our starting point is to 
linearize the original MB equations (1) around the SALT 
solution and to assess whether it is stable against small 
perturbations. In what follows we concentrate on the 
single-mode solutions only and thus insert the following 
expressions into the MB equations (1), 

E(x,t) = (E^x) +5E(x,t))e~ iuJl \ 

P{x,t) = (Pi(x) + 5P(x,t))e~ luJl \ (9) 

D(x , t) = (D(x) + 5D(x, £)). 

In this ansatz, Ei, Pi, D denote, respectively, the elec¬ 
tric field, the polarization and the inversion of the single¬ 
mode SALT solution of Eq. (4) and SE, 5P , 5D are the 
corresponding small perturbations around it. Utilizing 
the fact that the SALT solution also exactly solves the 
MB equations and neglecting the higher-order contribu¬ 
tions of the perturbations, we derive a set of linear PDEs 
(A.13) with respect to SE , SP, 5D. 

As a next step we convert the resulting system 
of equations into a standard eigenvalue problem. A 
split of the complex variables into their imaginary 
and real parts gives rise to a set of linear equa¬ 
tions for five independent fields, which for convenience 
can be represented as a single vector field, F(x) = 
(Re(AP), Im(AP), Re(AP), Im(AP), 5D). We look for so¬ 
lutions of the form F(x,t) = F(x)e at , where a is the 
growth rate, and derive a set of linear equations (A. 16) 
containing the spatial dependence only. Using an appro¬ 
priate discretization scheme and taking into account the 
periodic boundary conditions, we finally end up with a 
quadratic eigenvalue problem of the following form: 

AF + aBF + g 2 CF = 0, (10) 

where A, B and C are the corresponding matrices whose 
dimensions depend on the chosen spatial discretization. 

This eigenvalue problem can be solved numerically, 
resulting in a set of eigenvalues and eigenvectors, 
{a- 7 , Fi(x)}. Note that eigenvalues with Re(cr J ) > 0 
stand for the perturbations, which grow exponentially 
in time implying that our SALT solution is unstable. 
Conversely, if all eigenvalues have a real part smaller 
than zero, the SALT solution is stable against small per¬ 
turbations. Therefore, finding the eigenvalue with the 
largest real part is sufficient to classify the stability of 
the SALT solution. The imaginary part of a stands for 
the frequency relative to coq, with which the perturba¬ 
tion oscillates. In Fig. 2, we show a typical example of 
the eigenvalue spectra {cr- 7 } with different values of yy 
for the case of the symmetric ring laser described in the 
previous section. Note that due to the fact that the MB 
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equations for the single-mode regime are invariant under 
a global phase rotation, the value a = 0 always shows up. 
It, however, does not affect the behavior of the system 
and therefore is always excluded from the consideration. 

To assess the stability diagram of a SALT solution 
we start at a certain value of the pump strength and 
then gradually vary the value of yy. The dashed lines 
in the stability diagram depicted in Fig. lc correspond 
to the stability thresholds at which the eigenvalue with 
the largest real part, maxj [Re(cr J )], crosses the imagi¬ 
nary axis [Re(cr) = 0]. We emphasize that the bound¬ 
aries between stable and unstable regions, which we find 
in this way, are in excellent agreement with the time- 
dependent simulations, as seen in Fig. lc. It should be 
noted that previous studies on such a linear stability anal¬ 
ysis involved more restrictive approximations and a lim¬ 
ited class of perturbations by keeping track of low order 
Fourier terms only [32] , whereas our approach is exact in 
the framework of the MB equations. 

Note, that for the above procedure to work it is not 
necessary to compute the whole eigenvalue spectrum 
{cr- 7 }, but it is sufficient to only consider eigenvalues in 
the complex region close to the real axis and for imag¬ 
inary parts in the range of 0 to uj\. In systems where 
the rotating wave approximation is justified, only per¬ 
turbations within this frequency range can realistically 
influence the system. In systems with a complex refrac¬ 
tive index, spurious solutions that need to be excluded 
from the analysis, can occur in the region |Im(<r)| > uq. 
This restriction on the eigenvalues we are looking for can 
be used together with an iterative eigenvalue solver to 
check relatively quickly if a solution is stable or not [52] . 

It is worth noting that for the above stability analy¬ 
sis of the MB equations no additional assumptions have 
been made and the results are therefore valid for arbi¬ 
trary single-mode SALT solutions. In the next section 
we consider the ring laser with an embedded scatterer 
where modes are only near-degenerate to see how our 
results change for the case of such a slight lift of the de¬ 
generacy. 

EXAMPLE 2: ID RING LASER WITH BROKEN 
SYMMETRY 

In any real-world implementation of a ring laser or mi¬ 
crodisk laser, the rotational symmetry will be slightly 
broken due to inhomogeneities in the material or by im¬ 
perfections of the manufacturing process [30, 53]. While 
intuitively one would expect that a slight modification of 
the system should not change the structure of the lasing 
solutions, there is evidence that this symmetry breaking 
can have a strong effect on the lasing modes [53, 54]. In 
order to better understand the impact of such a slight 
symmetry breaking on the stability of the SALT solu¬ 
tions, we analyze here a ring laser where the rotational 



Figure 2. (Color online) Stability eigenvalues cr J of the ring 
laser for D 0 = 0.06 and yy = 10 -3 (blue circles), yy w7-10 -3 
(green crosses) and yy = 10 _1 (red Ys). All other parameters 
are chosen as in Fig. 1. One can see that the single-mode 
lasing solution is only stable for the intermediary value of 
7 y (green crosses), since for both other cases, a a j with a 
real part larger than zero exists. The stability eigenvalues 
with Im(cd) ~ 10 -2 are related to the competition between 
the nearly degenerate modes. The other visible clusters of 
stability eigenvalues at Re(cd) ~ —yy are related to damped 
relaxation oscillations observed when initially turning on the 
laser. The eigenvalues with Im(cr J )> 10 are related to other 
resonances of the system that could potentially also turn into 
active lasing modes. 


symmetry is broken by a scatterer as depicted in Fig. 3a. 

Since the inhomogeneity in this system breaks the ro¬ 
tational symmetry, solving Eq. (6) for the inactive system 
will not produce any traveling wave solutions. Instead, 
standing wave solutions similar to those found in the 
completely symmetric system are found. But while the 
standing wave solutions of the symmetric system always 
occurred in degenerate pairs, breaking the symmetry lifts 
the degeneracy such that the new modes slightly differ in 
their complex frequencies uJi. Therefore, they will nei¬ 
ther possess exactly the same lasing frequency, nor will 
they reach the lasing threshold at exactly the same same 
pump strength Dq. Since the frequency splitting does 
not fulfill the stationary inversion approximation for re¬ 
alistic values of 7 ||, a possible two-mode solution can not 
be described by SALT (at least not without an explicit 
stability analysis). 

In our analysis we consider gain parameters that only 
support the lasing action of a single pair of such modes 
(here located at UJi « 62) such that only the single-mode 
solutions of SALT seem viable candidates for steady-state 
lasing. These single-mode SALT solutions are shown in 
Fig. 3b in form of their frequency dependence on the 
pump strength Dq. In particular, the two single-mode 
solutions of SALT corresponding to the pair of nearly 
degenerate resonator modes are shown as mode A and 
B. When tracking mode A from its threshold while in¬ 
creasing the pump strength, it shows the following be- 
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Figure 3. (Color online) Results for a ring laser system as in 
Fig. 1 with 7 || =0.01, but with an additional scatterer added 
with l = 0.05 and yT = 1.05 + 0.0002L (a) Sketch of the sys¬ 
tem. (b) Frequencies of single-mode SALT solutions. It can 
be seen that the two degenerate modes of the symmetric sys¬ 
tem from Fig. 1 have split into modes A and B. Additionally, 
at Do ~ 0.003, two predominantly traveling wave solutions 
Ci and C 2 branch off from solution A. These two solutions 
are mirror-symmetric images of each other. Dotted lines mark 
unstable, solid lines mark stable solutions, (c) Intensity dis¬ 
tribution of modes A and C\ with the scatterer marked in 
gray. Mode A is shown directly at threshold (Do ~ 0.0015), 
mode Ci at Do = 0.1. It can be seen that while A is a stand¬ 
ing wave with minimal intensity 0 at the nodes, the predom¬ 
inantly traveling wave C 1 features a non-vanishing intensity 
everywhere in the cavity. The stability of mode Ci ,2 is studied 
in Fig. 4. 


havior: At D 0 ~ 0.0015 the standing wave single-mode 
solution A becomes active and the system starts lasing 
(see Fig. 3b). At this pump value the mode is stable, both 
according to the traditional SALT criterion for stability 
as well as in the linear stability calculations. However, a 
very small increase of the pump strength brings the eigen¬ 
value uji of the second mode of this near-degenerate pair 
to the real axis and thus renders the single-mode solution 
unstable (again, according to both criteria). Whether a 
stable two-mode solution in this near-degenerate regime 
exists is a question that presently falls outside the scope 
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Figure 4. (Color online) Stability of the SALT solution for 
the broken symmetry ring laser shown in Fig. 3. The ana¬ 
lyzed solution is the one marked as branch Ci ,2 in Fig. 3b. 
While the shape of the stability region is very similar to the 
one observed in Fig. 1, the broken symmetry system requires 
a minimal pump strength Do to compensate for the frequency 
splitting of the nearly degenerate modes before a stable pre¬ 
dominantly traveling wave solution can emerge. 
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of SALT (due to the non-stationary inversion) [39] . How¬ 
ever, with the techniques presented above we can investi¬ 
gate the stability of a single-mode solution in the regime 
beyond the point where the second mode passes the in¬ 
stability threshold. For this purpose, we track each of the 
two modes A and B towards higher pump strength using 
the SALT equation Eq. (4) for each mode as if the other 
mode was not active. Using, in addition, the stability 
criterion from the previous section reveals that both so¬ 
lutions on their own remain unstable for higher values of 
the pump strength. We find, however, that a bifurcation 
occurs at Do ~ 0.003, at which two further single-mode 
SALT solutions C\ ,2 of Eq. (4) branch off from solution 
A. These new modes share the same frequency, which 
lies approximately in between the frequencies of the un¬ 
stable modes A and B. To find these two modes one 
can, e.g., use a linear superposition of mode A and B 
at higher pump strengths as an initial guess for solving 
Eq. (4), similarly to the way the traveling wave solutions 
of the symmetric ring laser can be expressed as linear 
superpositions of the standing wave solutions. 

By looking at the mode profiles of solutions Ci ,2 {C\ 
is shown in Fig. 3c), we observe that these modes are 
related to the traveling wave solutions observed in the 
system with unbroken rotational symmetry. There, the 
two stable solutions were purely clockwise or counter¬ 
clockwise traveling waves. Here, each of the modes C\^ 
still features a dominant contribution in one direction 
and is identical to the other mode when being reflected 
at the symmetry axis containing the scatterer (x = 0.5). 
The major difference to the solutions observed in the un¬ 
broken ring laser is that the solutions Ci ,2 do not exist 
below a critical pump strength D cr i t ~ 0.003, since the 
nonlinear term in Eq. (4) needs to be strong enough to 
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compensate for the frequency splitting between the two 
modes. The fact that the nonlinear term is responsible 
for compensating the frequency difference between the 
two nearly degenerate modes was further investigated. 
We used a linear approximation to estimate the frequency 
shift experienced by the passive mode of the system that 
corresponds to the single-mode solution B as a conse¬ 
quence of the spatial hole burning from the active mode 
A. For this purpose, we modeled the change in the popu¬ 
lation inversion A D(x) as well as the resulting change in 
frequency A ujb of the passive mode as a linear function of 
the pump strength D 0 . The estimate predicted the two 
modes to share their frequency at a value of Do « 0.30, 
which is very close to the point Do ~ 0.33, where the so¬ 
lutions Ci ? 2 really emerge. The remaining small discrep¬ 
ancy can be attributed to the inaccuracy of the employed 
approximations. 

Next, we check the stability of the mode pair Ci ? 2 - 
First of all we recall that using the traditional SALT 
criterion Eq. (8) (which is not applicable here), one would 
find that the solutions Ci ? 2 are never stable. However, 
both the results of the FDTD calculations as well as the 
linear stability calculations do show that these modes are 
stable in well-defined limits. These limits are indicated in 
Fig. 4, where the stability of the single-mode solutions are 
depicted as obtained from both methods under variation 
of the pump strength Do as well as of the relaxation rate 
of the inversion yy. As discussed in the above paragraph, 
there is a small region of stability for very low values 
of pump strength Do where mode A is stable. Of more 
interest, however, is the large region of stable lasing for 
modes C\^ which opens up at the critical pump strength 
Dcrit ~ 0.003 (green region in Fig. 4) and which is similar 
to the one observed for the fully symmetric ring laser 
(Fig. 1). 

Altogether, the results of the linear stability analysis 
are in excellent agreement with the findings of the time- 
dependent simulations of the MB equations, as can be 
seen both in Fig. lb and in Fig. 4. The linear stabil¬ 
ity analysis thus proves to be a reliable tool to classify 
the stability of SALT solutions and can therefore be used 
even for systems where verifying the results through time- 
dependent simulations is not feasible. This is, in par¬ 
ticular, the case for higher dimensional systems, such as 
the microdisk system analyzed in the next section, where 
time-dependent simulations for a large range of parame¬ 
ters is computationally very demanding. 


EXAMPLE 3: 2D MICRODISK LASER WITH 
WEDGE 

Next, we will extend our results to a two-dimensional 
system, where full time-dependent simulations become 
impractical for the required range of parameters, but the 
linear stability analysis can still be easily performed to 



D 0 


Figure 5. (Color online) Single-mode lasing states for a 2D 
cavity with a wedge. The disk has a radius of 1, a refrac¬ 
tive index of n = 2 + 0.0 lz and the wedge on the right side 
of the disk has a depth of 0.05 and a width of 0.05. The 
gain parameters are uo a — 4.83 and yj_ = 0.1. (a) The cav¬ 
ity is outlined in white with arrows marking the position of 
the wedge. The spatial intensity pattern of two stable single¬ 
mode solutions are shown: (left) Standing wave mode B at 
pump strength Do = 0.0773 with regularly spaced nodes with 
intensity zero; (right) dominantly traveling wave solution C\ 
at pump strength Do = 0.15. (b) Laser frequencies of single¬ 
mode solutions to the SALT Eq. (4). Curves A (orange) and B 
(red) correspond to two standing wave solutions, which have 
even and odd symmetry with respect to the x-axis. The two 
modes represented by curve Ci ,2 (blue) feature a broken sym¬ 
metry, but are mirror-symmetric to each other. Solid (dotted) 
lines denote a stable (unstable) solution for qy = 0.01. Note, 
that the pump axis is separated into three regions of different 
linear scaling for the sake of clarity, (c) Stability diagram of 
the single-mode solutions shown in (b) with respect to the 
pump strength Do and the relaxation rate of the inversion, 
7 ||. The color-coding is as in Fig. lc except that the stability 
is here solely determined by the linear stability analysis. 


justify the stability of the SALT solutions. The system 
we will consider is a perturbed 2D microdisk laser where 



















the cavity is slightly deformed by cutting a small wedge 
into the disk (see arrow and white outline in Fig. 5a). In 
analogy to the scatterer for the ID ring laser, the wedge 
has the effect that the two originally degenerate threshold 
modes of the microdisk split (here the frequency splitting 
Auj ~ 1.4-10 -3 ). Hence, the multi-mode condition yy <C 
Auj cannot be satisfied for reasonable values of yy such 
that a two mode solution does not fulfill the SIA and the 
traditional SALT algorithm cannot be applied. 

The strategy to find the single-mode SALT solutions is 
a bit more sophisticated as compared to the ID system. 
We start with the threshold modes A and B correspond¬ 
ing to the nearly degenerate resonances of the unpumped 
system and track them while gradually increasing the 
pump strength Dq (see Fig. 5b). Mode B has a lower 
threshold and is stable in a tiny pump region starting 
from its threshold at D 0 ~ 0.0771 up to about 0.0777. 
The corresponding stable parameter region for mode F>, 
shown in the stability diagram in Fig. 5c by the green 
region on the very left of the figure, depends on both 
the pump strength Dq and on yy. After mode B has 
become unstable the mode is tracked further (neglecting 
the presence of mode A whose resonance eigenvalue has 
meanwhile also crossed the real axis). At Dq « 0.124 a 
pair of solutions C i ? 2 branches off from mode B. Whereas 
mode B and mode A feature a perfect even and odd sym¬ 
metry with respect to the x-axis (i.e., the symmetry axis 
of the system), the solutions C\^ do not possess this 
symmetry, but rather are mirror images of each other 
(compare mode profiles of modes B and C\ in Fig. 5a). 
In fact, the symmetry of the system is spontaneously bro¬ 
ken when either of these modes are lasing, a phenomenon 
that has been previously observed in simulations as well 
as in experiments [5]. 

In contrast to mode D, mode A is never a stable laser 
mode since it has a higher lasing threshold. At Dq « 
0.8 both (unstable) modes A and B feature two further 
branches. These can be understood as follows: Since 
the wedge of the two dimensional cavity only represents 
a small perturbation to the system, the symmetry with 
respect to the y -axis is only slightly broken, and, hence, 
modes A and B are nearly symmetric with respect to this 
axis (see left panel of Fig. 5a for the intensity pattern 
of mode B). At Dq ~ 0.8, this near-symmetry is no 
longer realized. However, since the modes have never 
been fully symmetric, there is not a single point at which 
the symmetry breaks, but rather a smooth transition (as 
compared, e.g., to the symmetry breaking transition with 
respect to the x-axis at Dq « 0.124). For mode L>, this 
smooth transition is clearly visible in Fig. 5b. For mode 
A this transition occurs in a much smaller pump interval, 
since mode A has a node directly located at the wedge 
of the 2D cavity and its symmetry is therefore only very 
slightly distorted. 

To obtain modes Ci 5 2 , we need to track them back¬ 
wards in pump strength starting at the branching point 


at D 0 ~ 0.124 (see Fig. 5b). Reducing the pump strength 
further, we observe that a sharp bend occurs in the fre¬ 
quency dependence of these two modes at Do « 0.86, 
where each of the modes Ci : 2 has evolved into a domi¬ 
nantly traveling wave mode (see right panel of Fig. 5a), 
similar to the nearly degenerate ID ring laser. Note that 
the small contribution traveling in clockwise direction 
can be observed as a modulation in the intensity pattern. 
Beyond this turning point, modes C\^ become stable in 
a large region of parameters Dq and yy as depicted by the 
central green area in the stability diagram in Fig. 5c. The 
easiest way to find the branch C \,2 numerically is to sum 
the fields of modes A and B with an additional relative 
phase of 7r/2 at a pump strength D o > 0.86 and use this 
as a guess for the nonlinear solver to converge towards 
one of the dominantly traveling wave modes. The solu¬ 
tion can then be tracked to uncover the whole branch 
of modes C\^ Using similar superpositions of already 
known modes as starting point for the nonlinear solver, 
we found several more branches of non-linearly induced 
single-mode SALT solutions within the plotted frequency 
range, albeit none of them are stable (these branches are 
not shown in Fig. 5c). 

The traveling wave solutions C\^ only remain stable 
until a pump strength of D 0 ~ 1.3. Here, the single-mode 
SALT solution becomes unstable due to the fact that an 
additional mode would start to lase. We also find that 
for values of yy < 2 • 10 -5 the modes C\^ are never sta¬ 
ble, which highlights again how important it is to take 
into account the value of yy for assessing the stability of 
a SALT mode. Traditionally in SALT single-mode lasing 
solutions were implicitly considered as stable (without 
considering yy) [17, 22], but these single-mode solutions 
were always only identified for the case where only a sin¬ 
gle one of the eigenvalues ay of Eq. (7) has a non-negative 
imaginary part. Our results show that single-mode SALT 
solutions can also exist for the case of multiple eigenval¬ 
ues uJi featuring a non-negative imaginary part. In this 
new situation, a stability analysis is, however, indispens¬ 
able for correctly interpreting the solutions of the SALT 
equation. When the results of the stability analysis are 
taken into account, the SALT equation allows us to accu¬ 
rately describe the steady-state of these systems without 
requiring any time-dependent simulations. 

CONCLUSION 

In this work we demonstrate that SALT can be used 
to describe the single-mode lasing regime of resonators 
with degenerate or near-degenerate mode pairs. Our ap¬ 
proach builds on a careful tracking of SALT modes in the 
non-linear lasing regime together with a linear stability 
analysis to judge the validity of the resulting solutions. 
The accuracy of the stability analysis itself was tested by 
a comparison with full time-dependent simulations based 
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on the Maxwell-Bloch equations, which shows excellent 
agreement in all cases. 

Our approach is ideally suited to treat microdisk 
whispering-gallery-mode lasers, which were previously 
difficult to simulate with SALT and often only accessible 
through time-dependent simulations or through strongly 
simplified models. Generally speaking, our work paves 
the way to study interesting non-linear phenomena, such 
as bifurcating solutions etc. within the efficient frame¬ 
work of SALT. 

One obvious direction for further study is the general¬ 
ization of our stability analysis to systems in the multi- 
mode lasing regime, which will be the aim of a subsequent 
paper [39]. 
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Appendix A: Derivation of the linear stability 
analysis 


We start from a solution {Pi (a;), ooi} of the SALT 
Eq. (4) at a given pump strength P 0 - From this, one 
can construct the polarization, Pi, and inversion, P, that 
show up in the MB Eqs. (4) via 


D{X) 1 + \TiEi(x)\ 2 
P\(x) = T 1 D(x)E 1 (x). 


(A.ll) 

(A.12) 


Using these quantities, we insert the ansatz (9) into the 
MB Eqs. 1. Using the fact that (Eie~ luJlt , Pie~ lu)lt , D) 
is a solution of the MB equations, we linearize the equa¬ 
tions with respect to the perturbations. This results in 
the following partial differential equations for the per¬ 
turbations of the electric field, the polarization, and the 
inversion 


eSE = X7 2 SE-SP + ojj(SP + eSE) + 2ioj 1 (5P + eSE) (A.13) 

SP = (i(coq — oj a ) — 7_l) SP — i , y±(E5D + SED) (A.14) 

SD = — 7|| SD + ^ (SE P* + Pi SP* - SE*P 1 - P* SP ). (A.15) 


While this system of equations is linear in the perturba¬ 
tion, it includes the terms SP * and £P*, which can not 
be expressed as linear combination of SP,SE,SD. In or¬ 
der to produce a completely linear system of equations, 
we therefore split the two complex fields and correspond¬ 
ing perturbations, as well as a possibly complex dielectric 
function e into their respective real and imaginary parts 
and consequently split Eqs. (A.13) and (A.14) into four 
real equations. This yields a linear system of five equa¬ 


tions with purely real terms. These correspond to the 
five independent fields Re(5P), Im(5P), R e(SP), Im (SP), 
SD , which for convenience we can summarize as a single 
vector field P. 

In order to analyze if a solution Pi of the SALT equa¬ 
tion (4) is stable we need to check if a perturbation exists 
which doesn’t relax back to the stable solution. Hence, 
we make an ansatz of the form F(x,t) = F(x)e crt which 
yields 


(SP r — eiSEi c r SE r ^cr 2 T 2 (SPi H - CiSE r V c r SEpuj\ <j — (SP r — CiSEi H - € r SE r ^j uj 2 — W 2 SE r — 0 
(8 Pi + CiSE r + e r SEi)cr 2 — 2 (SP r — CiSEi + e r 5E r )uicr — (SPi + e^P r + e r 8Ei)uj\ — \7 2 SEi = 0 

SP r a + 7_l (SP r — DSEi — SDE]) + (uq — uJ a )SPi = 0 
SPi<7 + 7_L(£Pi + DSE r + SDEl) — (uq — uj a )SP r = 0 
SDa + 7|| (SD + P*6Ei + SP r E} - SE r Pl - E^SPi) = 0, 


(A.16) 


where we have abbreviated Re(-) and Im(-) through the impose boundary conditions for the perturbation SE. For 
sub-indices r and z, respectively. In addition we need to the periodic ID ring laser we can simply assume periodic 
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boundary conditions, i.e., 5E(pc i e f t ) = AE(xright) and for 
the 2D system \75E r ^=°° i(oq — icr)SE e r . In a next step 
Eqs. (A. 16) are discretized using a suitable discretization 
scheme [22] which leads to a quadratic eigenvalue prob¬ 
lem that can easily be linearized. In our calculations we 
have chosen the finite element framework FEniCS [55] 
for discretizing Eqs. (A. 16) and used a perfectly matched 
layer for imposing the outgoing boundary conditions in 
the 2D case. For solving the linearized quadratic eigen¬ 
value problem we have used SLEPc [56]. 
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